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Ejection  of  Particles  From  a  Charged  Infinite  Strip 


1.  INTRODUCTION 

The  problem  of  the  emission  of  particles  from  separated  electrified  half  planes 
has  been  considered  by  Shima,  Kalman,  and  Carini.  1  The  electrode  geometry  they 
treated  is  perhaps  the  simplest  model  one  can  use  to  approximate  the  emission  of 
charged  particles  from  an  electrode  at  a  given  potential,  towards  a  second  electrode 
at  another  potential,  fixed  with  respect  to  the  first. 

In  this  report  we  consider  the  next  logical  step  in  modeling,  namely  an  infinite 
strip,  of  width  be  (Figure  1)  imbedded  between  two  infinite  half  planes,  ab,  cd, 
placed  symmetrically  about  the  origin.  We  set  the  width  be  equal  to  two,  and  the 
potential  difference  between  interior  and  exterior  electrodes  equal  to  unity.  As  will 
subsequently  be  shown,  there  is  a  difference  in  asymptotic  properties  between  the 
infinite  half  planes  treated  in  the  earlier  model  and  the  infinite  strip  treated  here. 
Because  of  this,  the  behavior  of  particles  ejected  from  the  plane  of  the  electrodes 
differs  not  only  quantitively,  but  qualitatively  as  well,  in  these  two  cases.  In  par¬ 
ticular,  for  the  half  plane  case,  every  particle  ejected  from  one  electrode,  with  the 
electric  field  in  the  positive  y  direction  above  it,  which  crosses  to  the  space  over 
the  other  electrode  with  the  electric  field  directed  in  the  negative  y  direction, 

(Received  for  publication  29  January  1982) 

1.  Shima,  Y.,  Kalman,  G. ,  and  Carini,  P.  (1981)  Particle  Trajectories  In  A 
Model  Electric  Field  I,  AFGL-TR-81-0167,  ATTA1W507. 


ultimately  returns  to  the  plane  of  ejection.  For  the  infinite  strip  case,  particles 
ejected  from  the  center  of  the  strip  generally  escape  and  do  not  return  to  the  ejec¬ 
tion  plane.  This  result  seems  to  be  largely  independent  of  both  angle  and  energy 
of  ejection.  As  the  point  of  emission  approaches  the  edge  of  the  strip,  particles 
are  returned  to  the  ejection  plane  for  a  wider  and  wider  set  of  initial  angles  and 
energies;  and  in  the  limit  of  an  infinitesimal  separation  between  the  point  of  ejec¬ 
tion  and  the  strip  edge,  the  results  of  Shima,  Kalman,  and  Carini  are,  as  would 
be  expected,  recovered. 


Figure  1.  Geometry  and  Coordinate  System  for 
Emission  From  a  Strip 


2.  THE  EQUATIONS 

The  potential  above  an  infinite  strip  of  width  2,  and  with  unit  potential  differ- 

2 

ence  between  strip  and  surrounding  plane  may  be  obtained  from  the  conformal 
t  ransf o  rmation 


w  s  in 


by  taking  the  imaginary  component  of  w  and  is  found  to  be 


V  «  Vj  +  1/ir  tan 


(1) 


2.  Churchill,  R.  V.  (1948)  Complex  Variables  and  Applications,  McGraw-Hill, 
New  York,  Appendix  11,  figure  19. 
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Here.  0  £  tan”1  £  tr .  Vj.  is  the  potential  of  the  plane  outside  the  strip,  and  Vj  +  1 
the  potential  inside,  on  the  strip. 

Taking  derivatives,  we  have  for  the  equations  of  motion 


(2) 

(3) 


or  in  non-dimensionalized  form. 


(4) 


(5) 


Here  m  and  e  are  the  particle  mass  and  charge  respectively,  the  dot  is  a  time 
derivative,  and 


If  we  set 


x  =  1  +  $ 


y  =  C 

where  £  and  £  are  small  quantities,  we  approach  the  (positive)  strip  edge,  and 
Eqs.  (4)  and  (5)  become 


Sr  ■  »/'  ttV  «6> 

dr  i  +  C 


A-  — T 

dr 2  r+r 


(7) 


which  aside  from  the  factor  1  h  and  reversals  in  signs  are  the  equations  of  Shima, 
Kalman,  and  Carlnl. 
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3.  TRAJECTORIES  AND  SEPARATION  CURVES 


We  now  solve  Eqs.  (4)  and  (5)  numerically,  using  a  Haming  method 
initial  conditions 


3 


for  the 


T 


0;  x 


V  y  =  °* 


dx 
3t  = 


2  UQ  cos  a. 


dy  . 

3t  " 


2  U  sin  a 
o 


(8) 


where  U  is  the  initial  particle  energy  and  a  the  angle  of  emission  (see  Figure  1). 
We  do  this  for  four  initial  coordinates. 


xQ=  0.0  (1  -  xQ=  1.)  yo=  0 

x  =  0.68  (1  -  x  2  0.3),  y„  =  ° 

o  o  o 

xQ=0.  90  (1  -  xQ=  0.  1).  yQ  =  ° 

xq=  0.968  (l-xo=0.03).  yQ=  ° 

for  three  initial  energies, 

U  =  0.01.  U  =  1.0,  U  =  100; 

o  o  o 

and  for  a  variety  of  initial  angles  a . 


The  results  are  shown  in  Figures  2,  3,  4,  and  5.  It  is  seen  that  for  the  most 
part  the  graphs  terminate  in  straight  lines  with  slopes  of  45°  on  the  log-log  plot, 
indicating  asymptotic  motion  in  a  force  free  field  and  an  ultimate  escape  from  the 
emission  plane.  It  is  only  for  Uq  =  1,  for  lower  emission  angles,  and  for  emission 
from  near  the  strip  edge  that  particles  are  returned  to  the  plane  of  the  strip. 

-5 

In  Figure  6  we  see  several  trajectories  for  UQ  =  1  and  1  -  x  =  1. 0  X  10  . 

Since  not  only  1  -  x  but  also  the  maximum  value  of  1  -  x  and  the  maximum  value 
o 

of  y  are  much  less  than  unity,  we  expect  this  case  to  revert  to  that  studied  by  Shima, 
Kalman,  and  Carini,  and  the  trajectories  in  the  Figures  are  quite  similar  in  form 
to  those  that  they  obtain. 

Although  the  graphs  shown  in  Figures  2  through  5  are  useful  as  illustrations 
of  the  forms  of  typical  trajectories,  they  contain  more  information  on  the  specific 
trajectories  shown  than  is  likely  to  be  of  general  interest.  As  a  means  of  presenting 


3.  Carnahan,  B. ,  Luther,  H.  A. ,  and  Wilkes,  J.  O.  (1969)  Applied  Numerical 
Methods,  John  Wiley,  New  York,  p  393  ff. 
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essential  information  on  a  wider  variety  of  initial  conditions,  graphs  of  complete 
trajectories  are  highly  inefficient.  A  more  condensed  form  of  information  display 
would  be  a  series  of  graphs  giving  the  x- coordinate  of  the  point  of  impact  (for  those 
particles  which  do  return)  as  a  function  of  the  three  initial  parameters,  xq,  Uq,  a. 
Although  such  graphs  are  easy  to  obtain,  they  still  require  substantial  amounts  of 
computer  time,  and  we  do  not  present  any  graphs  of  this  intermediate  type  here. 

A  display  which  contains  even  less  information  on  specific  trajectories,  but  is  most 
compact  in  presentation  is  a  series  of  "separation  curves",  each  one  characterized 
by  a  single  initial  coordinate,  xq.  With  reference  to  that  initial  coordinate,  any 
point  (UQ,  a)  which  falls  above  the  designated  curve  represents  an  initial  condition 
corresponding  to  a  particle  which  ultimately  escapes  from  the  plane  of  ejection 
(y=  0).  Conversely,  any  point  falling  below  the  designated  curve  represents  an 
initial  condition  leading  to  a  trajectory  which  ultimately  returns  to  the  ejection 
plane.  From  this  single  figure  then  (Figure  7)  the  behavior  of  particles  with  regard 
to  ultimate  emission  or  return,  and  for  all  initial  conditions  xq,  U  ,  a  can  be 
inferred.  For  xq  <  0.7,  all  particles  escape  regardless  of  initial  energy  or  initial 
angle.  For  particles  0. 7  <  xq  <  0. 97  only  particles  with  low  angles  and  a  limited 
energy  range  UQ  falling  roughly  between  0.  1  and  -  10  are  returned.  As  xq 
approaches  the  edge  of  the  inner  strip  (xq  >  0.  97 )  particles  of  increasingly  higher 
initial  angles  and  broader  energy  ranges  are  returned  to  the  plane  of  origin. 
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Figure  2.  Typical  Trajectories  for  y0  =0,  xQ  =  0.0.  The  numbers  in 
parenthesis  are  emission  energy,  UQ,  and  emission  angle,  a 
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Figure  6.  Typical  Trajectories  for  Emission  From  Immediate  Vicinity  of  Strip 
Edge  Showing  Similarity  in  Form  to  Those  Obtained  by  Shima,  Kalman,  and 
Carini 


4.  THE  LIMIT  ALPHA  APPROACHES  ZERO 

It  is  clear  from  the  nature  of  the  problem  that  if  initial  conditions  yQ  >  0, 

<*  >  0  are  chosen  (that  is,  the  particle  is  fired  down  on  the  outside  electrode)  it  will 
always  be  possible  to  have  the  particle  reach  the  plane  of  the  electrode.  Hence,  it 
should  not  be  a  cause  for  surprise  that  solutions  show  a  unique  behavior  for  yQ  =  0 
and  a  ■*  0.  Although  values  of  a  plotted  in  Figure  7  are  only  accurate  to  within 
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±  0.  5  degree,  the  apparent  peculiarities  which  arise  for  small  angles  of  ejection 
were  investigated.  These  questions  are  more  of  theoretical  than  of  practical  inter¬ 
est,  since  from  an  experimental  viewpoint  there  is  no  possibility  of  discriminating 
between  ejection  angles  of  0  and  1  degree.  The  results  of  this  investigation  are 
presented  in  the  following  table.  Here,  the  numbers  given  are  the  x-coordinate 
for  which  the  particle  returns  to  y  =  0.  All  values  are  for  a  =  0.  0  degree. 


u 

o 

70 

80 

90 

100 

X 

o 

0.4 

8.  5 

8.  2 

8.  6 

8.  1 

0.5 

5.  8 

5.8 

5.6 

5.4 

0.6 

4.  0 

4.0 

4.0 

3.9 

0.7 

2.  86 

2.  8 

2.9 

2.8 

For  identical  values  of  Uq  and  xq,  but  with  a  =  1.  0  degree  the  particle  escapes 
for  all  Uq,  xq.  Thus,  there  is  a  radical  difference  in  behavior  between  initial 
conditions  having  a  =  0.  0  and  a  =  1.0. 


Figure  7.  Separation  Curves  for  Strip  Problem. 
For  any  value  (UQ,  a )  falling  below  a  given  curve, 
the  particle  is  returned  to  the  plane  of  emission; 
for  any  value  (U0.  a)  falling  above  a  given  curve, 
the  particle  escapes 
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5.  ASYMPTOTIC  PROPERTIES 

If,  in  Eqs.  (4)  and  (5)  we  set  p2  =  x2  +  y2  and  allow  x,  y  ■*  «  we  obtain 

lim  ^-4  ~  lim  ~  1/p2  .  W 

x,  y-»«  dT  x,  y-*»  dr  ^ 

2  2  2 

On  the  other  hand,  if  in  Eqs,  (6)  and  (7)  we  let  p  =  I  +  C  and  allow 
we  have 

lim  ^-1  ~  lim  -^-4  ~  1/p  .  dO) 

|,  C  -* «  dT  i.  ?  -* 00  dT  ( 

Thus  it  is  seen  immediately  that  the  two  sets  of  equations  have  different 
asymptotic  properties. 

To  study  the  effects  on  trajectories  of  the  asymptotic  properties  of  Eqs.  (4) 
and  (5)  we  seek  solutions  of  the  type 


x  =  or  +e 


y  =  0  T  +7} 


where 


‘1”!  sr  ■  lb  ■ 


OO 


7  oo 


If  we  put  Eqs.  (11)  and  (12)  into  Eqs.  (4)  and  (5),  we  have,  to  lowest  order 

a§ 


d2c 


1 


-  4/7T  - 2 - JT-B - 3- 

dr  "  (a  +0V  t  2 


d2q  _  a2  -  0  2  1 

— -4-  -  -  2/jt  — 2 - 9T  — r 

dr2  (a+02)2  t  2 


Imposing  the  boundary  conditions 

3?  =  ° 


(ID 

(12) 

(13) 


(14) 


(15) 
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we  find  that  £,  J|,  -  inr  and  hence  that  solutions  of  the  form  Eqs.  (11)  and  (12), 
approaching  the  limits  of  Eq.  (13)  are  possible.  This  is  consistent  with  the  results 
of  the  numerical  solutions. 

If  on  the  other  hand,  with  x,  y  replaced  by  |,  £  respectively,  we  put  Eqs.  (11) 
and  (12)  into  Eqs.  (6)  and  (7)  we  find  to  lowest  order 


d2e  = 

1 

(16) 

dr  ? 

T 

^4  = 

&r 

-a 

1 

(17) 

T  - 

Hence,  regardless  of  what  boundary  conditions  are  imposed  on  -  and  ^3  we  obtain 

lim  t ,  tj  ~  t  i  n  r  ,  and 
x  -*  00 

lim  =  lim  in  r  -*  «  . 

T  -*  «  a~  T  -*  «o 

Thus  there  cannot  exist  solutions  of  the  type  of  Eqs.  (11)  and  (12)  having  limits  of 
Eq.  (13).  This  is  consistent  with  the  results  of  Shima,  Kalman  and  Carini  who  find 
that  all  emitted  particles  return  to  the  plane  of  origin. 

6.  COMMENTS  ON  THE  CIRCULAR  EMBEDDED  ELECTRODE 

From  a  modeling  standpoint,  it  is  clear  that  a  circular  electrode  at  one  poten¬ 
tial  embedded  in  an  infinite  plane  at  a  second  potential  will  be  a  more  realistic 
representation  for  most  physical  electrode  systems  than  are  the  two  models  so  far 
discussed.  The  price  of  this  fidelity  is  increased  labor.  The  solution  of  the  bound¬ 
ary  value  problem  of  the  circular  electrode  is  available,  but  only  as  an  infinite 
integral, 

00 

V=V.+  J  J,(A)  e'Az  J  (Ap)  dA  .  (18) 

Here,  z,p  are  the  usual  cylindrical  coordinates,  JQ  and  J ^  are  zeroth  and  first  order 
Bessel  functions,  and  a  unit  potential  between  inner  and  outer  electrodes  is  again 
assumed.  In  the  plane  of  the  electrodes  (z  =  0)  we  have  the  infinite  integral 
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(19) 


00  1  (J^  <  1 

j  J  AX)  J  (p  A)  dA  =  1/2  p2  =  1  . 

0  0  ft*  >  1 

Along  the  axis  of  symmetry  the  integral  is  expressible  as  an  elementary  function 
00 

/  J,(A)  e'Xz  dA=  1 - 1  z  (20) 

•'o  1  yl+z7 

4 

The  general  integral  Eq.  (18)  has  been  tabulated  for  a  limited  range  of  values  of 
z  and  p.  To  compute  trajectories  however,  it  is  necessary  to  follow  particles  out 
to  large  values  of  the  independent  variables.  This  might  be  accomplished  by  using 
a  multipole  expansion  for  large  values;  either  a  numerical  solution  or  tabulation 
for  small  values;  and  matching  at  the  boundary  between,  all  in  all,  requiring  a  fair 
amount  of  labor. 

The  labor  is  probably  not  justified  because  of  the  following  considerations. 
From  Eq.  (20),  the  asymptotic  properties  of  the  disk  problem  are  available.  We 
have 


lim  V  ~  \ 
z  -*  <*>  tT 

and  for  the  field 

lim  E  ~  — j  « 

Thus  the  two  half  planes,  the  strip,  and  the  circular  disk  form  a  sequence  with 
respect  to  their  asymptotic  behavior.  We  expect  they  will  also  form  a  sequence 
with  respect  to  the  behavior  of  particles  emitted  from  their  surfaces.  Since  most 
initial  conditions  lead  to  particles  which  escape  from  the  strip,  this  will  be  even 
more  the  case  for  the  circular  disk.  We  expect  that  only  under  very  restricted 
initial  conditions  will  particles  emitted  from  a  circular  disk  return  to  their  plane 
of  origin.  The  labor  that  would  be  required  to  solve  for  trajectories  in  this  case 
is  thus  probably  not  justified. 
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